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ABSTRACT 

We have studied the star formation history and the initial mass function 
(IMF) using the age and mass derived from spectral energy distribution (SED) 
fitting and from color-magnitude diagrams. We also examined the physical and 
structural parameters of more than 1,000 pre-main sequence stars in NGC 2264 
using the on-line SED fitting tool (SED fitter) of Robitaille et al. 

We have compared the physical parameters of central stars from SED fitter 
and other methods. The temperature of the central star is, in many cases, much 
higher than that expected from its spectral type. The mass and age from SED 
fitter are not well matched with those from color-magnitude diagrams. We have 
made some suggestions to improve the accuracy of temperature estimates in SED 
fitter. 

In most cases these parameters of individual stars from SED fitter in a star 
forming region (SFR) or in the whole cluster showed nearly no systematic vari- 
ation with age or with any other relevant parameter. On the other hand, the 
median properties of stars in NGC 2264 SFRs showed an evident evolutionary 
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effect and were interrelated to each other. Such differences are caused by a larger 
age spread within an SFR than between them. 

The cumulative distribution of stellar ages showed a distinct difference among 
SFRs. A Kolmogorov-Smirnov test gave a very low probability of them being 
from the same population. The results indicate that star formation in NGC 
2264 started at the surface region (Halo and Field regions) about 6-7 Myr 
ago, propagated into the molecular cloud and finally triggered the recent star 
formation in the Spokes cluster. The kind of sequential star formation that 
started in the low-density surface region (Halo and Field regions) implies that 
star formation in NGC 2264 was triggered by an external source. 

The IMF of NGC 2264 was determined in two different ways. One method 
used the stellar mass from the SED fitting tool, the other used the stellar mass 
from color-magnitude diagrams. The first IMF showed a distinct peak at m « 
2M Q , but the second did not. We attribute the peak as an artifact of the SED 
fitting tool because there is no such IMF with a peak at m « 2M . The slope 
of the IMF of NGC 2264 for massive stars (logm > 0.5) is -1.7 ± 0.1, which is 
somewhat steeper than the so-called standard "Salpeter-Kroupa" IMF. We also 
present data for 79 young brown dwarf candidates. 

Subject headings: open clusters and associations: individual (NGC 2264) - stars: 
luminosity function, mass function 



1. INTRODUCTION 



The initial mass function (IMF) is one of the core topics in st ellar astronom ical research. 
The universality of the IMF has been debated for a long time (see IScalol (120051 ) for example). 
In many cases the slope of the IMF of open clusters is well consistent with the so-called 
standard "Salp eter-Kroupa" IM F within the observational uncertainties , eg. T = —1.2 ±0.1 



for NGC 6231 (ISung et al.lll998l ). Tr 14 an d Tr 16 in the Carina nebula flHur et al.ll2010l ). or 
the Arches cluster in the Galactic center (jEspinoza et al.l 120091 ) . But in some cases, a top- 
heavy IMF i s repo rted for massive young open clusters such as Westerlun d 1 (r = — 0.8±0.1 - 
Chun et al.l (120101 ) ) or for young stars in the Galactic center stellar disk (IBartko et al.1120 10'). 
On th e other hand, a deficit of massive stars has also been reported for Taurus (ILuhman et al. 



20031 ) and the Cone nebula region of NGC 2264 (r = -1.7 ± 0.3 - [Park et al.l ( 12000( 1). In 



chemical evolution models o f bulges of our Galaxy and M31, a non-standard IMF has been 
favored ( iBallero et al.l 120071 ): however, recent precise abundance analyses of bulge and disk 
stars favors a universal IMF (G. Gilmore - private communication). The universality of the 
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IMF remains an intriguing problem to be solved in the near future. Another outstanding 



issue in IMF research concern s the companion mass function in binary systems (jLadall2006 
Metchev & Hillenbrand! bopgh . 



The origin of the IMF is a hot topic of research. The similarity between the dense core 
clump mass function and the stellar IMF has led many investigators to implicitly assume 



that the clump m ass f unction is the origin of the stellar mass function (IMotto et al.lll998 



Alves et al.l 120071 ) . But ISmith et al.l ( 120091 ) show the non-existence of such a correlation be- 
tween clump mass and stellar mass. Owing to the rapid development of computing power 
and techniques many attempts have been made to understand the origin of the IMF through 
numeri cal modeling. N owadays vario us physical processes such as gravitati onal fragmen- 



tatio n (IZinnecker 



984 ), turbulence (IPadoan &: Nordlundl 120021 ) . accr etion (Bonnell et al. 



20011 ) . coalescence (jElmegreen fc Shadmehril 120031 ) , and mag n etic fi elds ( jShu et al.ll2004l ) are 



taken into account in the simulations. Recently iDib et al.l (120101 ) were able to reproduce 
IMF or clump mass functions nearly identical to observations. 

Star formation history remains one of the most interesting topics of current research. It 
is also related to triggering mechanisms of star formation. Most Milky Way OB associations 
extend over a few degrees (linear size of about 100 pc). It is very difficult to imagine that 
star formation on such a large scale occurred at a single instant. We can therefore expect to 
see age differences within an OB association. Classical examples of this are the sequential 
star formation in the OB associations - Sco-Cen associa tion (= Upper Scorpius + U pper 
Centa urus and Lup us + Lower Centaurus and Crux) (IPreibisch &: Zinneckei Il999[ ). Ori 



OBI (Blaauw 



Wright et a 



1964 ) or the W3-W4-W5 regions in Cassiopeia (ILada et al.lll978l ). Recently. 



(201 01) ide n tified two different age groups in the massive OB association Cyg 



OB2, while iBillot et al.l ( 120101 ) could not find any sig nature of propagating star formation 
in the Vulpecula OB association. On the other hand, iPovich fc Whitney! ( 120101 ) suggested 
a possible relation between star formation history and the IMF. They claimed evidence 
for delayed massive star formation in the M17 SWex region from deviations in the young 
stellar object (YSO) mass function from the Salpeter-Kroupa IMF. Another issue relating 
to the IMF and star formation history is the IMF of field stars such as the stars in the 
solar neighborhood. Without detailed knowledge of the star formation history in a region 
it is nearly impossible to trace the underlying IMF from the present day mass function (see 
Elmegreen fc Scald (120061 ) for example). 



This is the final paper based on extensive optical photometry of NG C 2264. The IMF 
aroun d S Monocerotis and for low- mass s tars in NGC 2264 was derived in ISung et al.l ( 120041 . 
2005h (Paper I & II). ISung et all (boost ) (Paper III) published an extensive photometric 
catalogue of NGC 2264. They identified two active star forming regions (SFRs - S MON & 
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CONE) and a less active surrounding halo region (HALO). More recently, ISung et al.l (120091 ) 
studied the spatial and spectral distributions of yo ung pre-main seq uence stars in NGC 2264 
based on Spitzer IRAC (Infrared A rray Camera - (Fazio et al.l 120041 )) and MIPS (Multiband 
Imaging Photometer for Spitzer - ( Rieke et al. 20041)) 24/xm observations. From the spatial 
distribution of Cl ass I and II YSOs, ISung et al.l ( 120091 ) divided the CONE nebula region of 
Sung et al.l (120081 ) into three SFRs - the Spokes cluster, Cone(C) and Cone(H). 



This paper is organized as follows. In §2 we describe input data and some results from 
the SED fitting model (hereafter SED fitter). The star formation history and the IMF of 
NGC 2264 are presented in §3 and §4, respectively. We estimated the age and mass of stars 
in NGC 2264 in two ways, one from SED fitter, the other from color-magnitude diagrams 
(CMDs). In §5 we present a catalogue of young brown dwarf candidates (YBDCs) and 
their selection criteria. T he effect of binarity on the IMF and the IMF of Ml 7 SWex by 
Povich fc Whitney! (120101 ) are discussed in §6. The summary and conclusions are in §7. 



SPECTRAL ENERGY DISTRIBUTION MODEL FITTING 



In order to characterize the physical properties and the evolutionary status of disks and 
envelopes around pre-main sequence (PMS) stars in NGC 226 4, we fitted the SEP of stars 
detected in the mid-IR Spitzer observations to the model grid of iRobitaille et al.l ( 120061 . 120071 ) 
(we call this web pag e the SED fitter. ) j] : To obtain better results we comp i led and used all 



the a vailable optical (ISung et al. 



19971: [Park et al.l I200CH : ISung et al.l 120081: iDahm & Simon 



20051 ). near-IR 2MASS flSkrutskie et al.l 120061 ). and mid-IR Spitzer data (ISung et al.ll2009f ). 
We limited our attention to objects which were detected in all IRAC bands and whose total 
photometric error in the four IRAC bands was smaller than 0.25 mag. The number of 
targets was 1,173 (381 stars in S Mon, 185 stars in the Spokes cluster, 157 stars i n Cone(C), 



206 st ars in Cone(H), 178 stars in Halo, and 66 stars in Field regions - see ISung et al. 



( 120091 ) for the region names). Among them, 853 objects were known members (X-ray, Ha 
emission, IR excess from Spitzer IRAC, or early type stars) of NGC 2264. In addition, 
33 embedded objects in S Mon, Spokes cluster, Cone(C) and Cone(H) were also included. 
Among the above 886 objects, 336 objects (38 %) were detected in MIPS 24 \im which 
gives an important constraint in SED fitting especially for the YSOs with transition disks or 
with pre-transition disks. The total number of known members in NGC 2264 is 1,162 stars, 
and therefore the fraction of stars used in SED fitting was 73.4 % (= 853 / 1162) of the 
total known members. In addition we also examined physical properties of 40 BMS (below 



1 http : / / caravan . astro . wise . edu/protostars/ sedf itter .pho 
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pre-main sequence) stars. 

SED fitter gives a wide range of models for a given SED. In such cases we should have 
model selection criteria. If several models give similar SEDs for a target, it is best to select 
a model that gives the smallest mass for the central star because (1) low-mass stars have a 
higher probability of having disks or envelopes, and (2) low-mass stars occur more frequently 
than massive stars. In addition, in many cases, SED fitter gives a higher temperature than 
the temperature derived from the spectral type of the star (see §2.1.1). In most cases the 
temperature of the Xmin m °del is equal to or higher than that of the mass-minimum model. 
The best model selection criterion is therefore, the model which gives the smallest mass for 
the central star among the 10 best models from SED fitter if x 2 1 Xmin — 2-0 or \ 2 < 2.0 for 
the case of Xmin — 1-0- But the best model for a BMS star is, on the other hand, the model 
whose x 2 is the minimum among the candidate models because if BMS stars are really PMS 
stars with edge-on disks, the criterion used for normal stars above may give a biased result. 

In SED fitting we applied constraints on the distance and reddening, logd (kpc) = -0.15 
-0.0 and Ay = - 100 mag. The parameters from SED fitter are those for the aperture of 
2762 AU (about 3."5 at 760 pc). 



2.1. Comparison of Stellar Parameters 

The physical or structural properties of PMS stars evolve as the central star ages. We 
therefore expect an evolution or variation of such properties. Young open clusters are ideal 
targets for this kind of study and it is important to check these properties with the SED fitting 
tool. The results can give some insights into an understanding of the structural evolution of 
PMS stars as well as the evolution of the cluster itself. In addition, it also provides evidence 
on the reliability of the SED fitting tool, and thereby contributes to improving SED fitter. 

In this section we compare several parameters relating to the central star from SED 
fitter with those from other methods. 



2.1.1. Effective Temperature 



Spectral type is the primary indicator of stellar temperature. iRebull et al.l ( 120021 ) pub- 
lished the spectral type of 360 stars in NGC 2264. Among them 174 stars are targets for SED 
fitting. Figure [T] shows the relation between spectral type and temperature. The left panel 
shows the relation between spectral type and the temperature from SED fitting, while the 
right panel shows the same relation but with temperature from the color versus temperature 
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relation of iBesselll ( 119951 ). As there is no known way to determine the reddening of late type 
PMS stars u sing photometric data alone, we app lied the same reddening E{B — V) = 0.07 
for all stars ( Sung et al.lll997t IPark et al.ll200C ) □■ The solid line in the figure is the spectral 
type versus temperature relation of MS stars ( ISchmidt-Kalerlll982[ ). 



Normally the expected error in spectral classification is 1 - 2 subclasses, but as many 
of our targets are Class II objects we can expect an error of 2 or more subclasses in spectral 
type due to the intrinsic variability In the left panel of Figure [TJ there are many stars with 
much higher temperatures from SED fitter than that from the spectral type. Although we 
take the error in spectral type into account, the difference in temperature is excessive. The 
reddening of these stars would be more than 2.4 mag (mean Ay = 3.32 mag). This value 
is much larger than the median Ay of all SFRs except Spokes cluster in Table [TJ This 
means that the reddening from SED fitter is less reliable, and may cause a problem in the 
temperature estimate. 

On the other hand the temperature derived from color indices is relatively well matched 
with that expected from its spectral type. There are three exceptions - C35587, C39168, and 
C36097D (= C36097 + C36109). Two of them are due to a large error in spectral type (G6: 
for C35587 and K0: for C39168). C36097D may be due to duplicity or reddening correction 
(SED fitter estimates the Ay of C36097D as 6.26 mag). 

The physical properties of the central source, such as effective temperature and lu- 
minosity, are the most important parameters in understanding the evolutionary status of 
the central source, and play a crucial role in the physical properties of disks and envelopes 
because dust in the disks and envelopes are heated by irradiation from the star and from 
accretion shocks at the stellar surface. A large difference in temperature of some stars in 
Figured] implies that SED fitter seems not to use optical color as a temperature indicator. As 
the optical color is a good temperature indicator, it is better for SED fitter to estimate the 
temperature of the central source using optical color, such as V — I which is least affected by 
UV excess from accretion shocks or IR excess from dust in the disks. In addition, color is, in 
most cases, obtained from semi-simultaneous observations, and therefore is less affected by 
the variability. As optical data, near-IR 2MASS data, and Spitzer mid-IR data are obtained 
at different epochs, there are many cases showing a systematic difference in fluxes among 
photometric bands due to variability. 

First determine the temperature of the central star from color indices, estimate the 
reddening Ay from SED fitting, and then estimate the temperature again using reddening- 



2 We also applied a different reddening for a different SFR (the median value of reddening from SED fitter 
in Table [T]), the result was very similar. 
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corrected optical colors. This iterative process seems to be required to get better results. 



2.1.2. Mass and Age 



Unlike stars in the main sequence (MS) stage, the mass, as well as the age of PMS stars 
can be determined from the CMD with the aid of theoretical evolution models of PMS stars. 
As there is no reliable method of estimating the reddening for low-mass PMS stars without 
prior knowled ge of spectra l type , we have applied th e mean reddening E(B — V) — 0.07 mag 
determined bv lSung et all f)l997h : IPark et al.l (boooh using early type MS stars in NGC 2264. 
This value is somewhat smaller than the m edian Ay from SE P fitter or E(B — V) from 
spectral type versus intrinsic color relation (iRebull et al.ll2002l ). But optically observable 
stars are relatively less embedded an d the difference i n age or mass resulting from a smaller 
value of reddening should be minor (jSung et al.l 120041 ) . 



The PMS e v olutio n models used are those by ISiess et al.l (120001 ) for m > 0.1 M^ and 



bvlBaraffe et al.l (119981 ) for m < O.1M (see Figure [9] of this paper or Figure 4 of ISung et al. 



(120041 ) ). We compared the mass and age from CMDs with those from SED fitter in Figure 
[2j The size of symbols is proportional to the amount of reddening Ay from SED fitter. 
Circles and squares represent Class II YSOs and stars with normal photosphere (or Class 
III), respectively. In the right panel open symbols denote massive stars (m* SED > 1 M ), 
while filled symbols stand for low-mass stars (m^sED < 1 M ). 

Masses of less reddened stars (small symbols) are well consistent with each other, but 
masses from SED fitter for high A^sed objects are in general larger than those from CMDs. 



And we easily see the lower mass limit 



m 



0.1 M Q ) of SED fitter which is the mass limit 



of PMS evolution models of ISiess et al.l (120001 ). In addition, there are two clumps in the 
distribution of masses from SED fitter - a distinct clump at logm w 0.2 - 0.3 and a less 
prominent at logm ~ -0.2 - -0.4. These two clumps are directly related to the peaks in the 
IMF of Figure 

The ages from SED fitter and from CMDs do not show any correlation (see the right 
panel of Figure [2]). For low mass stars (m < 1 M Q ), ages from SED fitter are in general 
smaller than those from CMDs, but for massive stars, ages from SED fitter are, on the other 
hand, larger than those from CMDs. The ages of massive stars from SED fitter are mostly 
larger than 5 Myr with an upper limit of 10 Myr, but those from CMDs are distributed 
widely from ~ to 13 Myr. For low-mass stars there is a strong concentration of ages from 
SED fitter between 0-3 Myr. But the ages from CMDs are distributed evenly between - 
6 Myr. Many Class II objects have relatively smaller ages from SED fitter. This fact implies 



-8 - 



that the age difference may be partly related with the reddening correction. But even when 
we applied the median value of reddening for a given SFR the correlation does not improve 
noticeably It is very difficult to surmise the source of such differences. 



2.2. SED Fitting Results for Individual Stars 

As we are dealing with many PMS stars in the young open cluster NGC 2264 and expect 
that many parameters relating to the structure or evolution of PMS stars are interrelated, it 
is therefore of interest to check whether there are any systematic variations of one parameter 
against another. 



2.2.1. General Properties 



The mean value of the distance (logrf, d in kpc) from SED fitting was -0.089 ± 0.051 (d 
= 815 ± 95 pc) for all PMS stars and the median value was between -0.09 and -0.11 (see Table 
[T]). This value is c onsistent with oth er determinations - 760 ± 90 pc from the zero age main 
sequence relation (ISung et al.lll997l ) and 913 ± 40 (sampling) ± 110 (systematic) pc from 
the assumption that the distribution of inclinatio n angles determin ed from the projected 
rotational velocity and rotation period is random ([Baxter et al.l 120091 ) . 



We could not find any correlation between age and reddening. This means that the 
interstellar material in NGC 2264 may be very inhomogeneous. The total luminosity did 
not show any systematic variation with age. This is a result of a combination of the mass 
spectrum and different evolutionary stages. If we limited our attention to the stars with 
masses less than 1 M Q , we could find a weak evolutionary trend between age and luminosity. 
The total luminosity of low-mass stars decreased with age with the power of -0.28 ± 0.02. 
Other parameters did not show any systematic variation with age or mass. 



2.2.2. Disks 

We present several figures relating to the disk properties in Figure [3j We show the 
relation between stellar age and disk accretion rate in the first panel of Figure |3J There 
seems to be virtually no variation in the upper limit of the disk accretion rate, but the lower 
limit decreases with time. The scatter in the disk accretion rate of \ogr age < —1 is about 3 
orders of magnitude, rising to about 8 orders of magnitude at logr affe ps +1. The relation 
between disk mass and stellar age shows a similar variation with time and a similar scatter. 



As a result, we see a rather tight correlation between disk mass and disk accretion rate. 



\og{M disk /M+) = -5.38(±0.06) + 1.09(±0.01) • \og(M disk /M+), (r = 0.936). 

There is an order of magnitude scatter around the relation at a given disk mass. In addition, 
the relation shows a slight curvature - higher accretion rates for massive disks. Stars in the 
younger SFRs - Spokes cluster and Cone(C) - are more concentrated in the upper right of 
the figure, i.e. higher accretion rates and more massive disks. 

The Rdisk,out of individual stars increases with age. The power of the lower boundary is 
about 0.47 and that of upper boundary is about 0.89. The initial scatter ranges over about 
an order of magnitude, but it increases with age and reaches about two orders of magnitude 
ne ar logr n , gP . ~ 1. And there are several outliers between logr a9e = - 1 as in Figure 6 



of iRobitaille et al.l ( 120061 ). Some of these outliers are double stars, which is reminiscent of 
disk truncation due to the existence of companion, but it is uncertain why they have smaller 
Rdisk,nut than the othe r YSOs. There is an exception which cannot be found in Figure 6 of 



Robitaille et al.l (120061 ) . C36015 has very small Rdisk,out- C36015 is a Class I object close to 



the brightest IR source in NGC 2264 IRS 1 (0/362 = 0.08 pc at 760 pc). 

If disk evolution occurred in a gradual manner, we could expect an increase in the 
disk inner radius with age. But the inner disk radius (Rdisk,in in units of R su b or AU) for 
individual stars shows no systematic variation with age. The disk scale height at 100AU, on 
the other hand, is greater for younger stars with massive disks (log z^xooau — 0.84(±0.02) + 
0.04(±0.00) -log^Mdisk/M*), r = 0.374). Other parameters, such as scale height factor Zf actor , 
flaring factor and disk accretion power a acc show no appreciable sy stematic variation with 



time. There is a very small variation of Zf act0 r and /3 with time, and IRobitaille et al.l (120061 ) 
interpreted such variations as the settling of dust to the disk mid-plane. 



2.2.3. Envelopes 

We present several figures relating to the envelope parameters in Figure HI The mass 
in the envelope is about 10 _1±1 M Q at the initial stage and decreases gradually with time 
reaching about 10~ 4±1 M Q at about 1 Myr. After that, the envelope masses show a large 
scatter. The mass in the envelope disperses rapidly and accretion from the envelope is halted 
at this stage (see the last figure of Figure H|). However, the large scatter in the envelope mass 
at r age > lMyr could simply be an artifact of the SED fitting. 

Although there is a large scatter, older stars have a larger cavity opening angle (log 4> cav ity 
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1.73(±0.01) + 0.33(±0.01) ■ logr a5e , r = 0.832). The cavity opening angle is about 10° at 
0.01 Myr, and increases up to about 55° at 1 Myr. The cavity density, on the other hand, 
decreases as age increases (log p ca vit y ~ — 20.93(±0.02) — 0.67(±0.02) logr a9e , r = —0.839). 
The initial value of the cavity density is about io~ 19 - 5±0 - 5 g cm -3 , and decreases with age to 
the power of -2/3, i.e. p C avit y oc r ag e^ 3 . At an age of 1 Myr, the cavity density approaches 
about io~ 210±0 - 5 g cm -3 . This fact implies that the sweeping up of envelope material occurs 
at a very early stage. 

The circumstellar extinction Ay fiircum (the extinction from the edge of the circum- 
stellar envelope to the surface of the central star) decreases with stellar age as Ay jCircum 
oc T~^ 75±om . There are three distinct groups in the diagram of envelope mass and cir- 
cumstellar extinction. The first is the highly extincted objects with massive envelopes, 
the second group is less extincted with less massive envelopes. The second group show 
a good correlation between envelope mass and circumstellar extinction (log Av jC ircum ~ 
-1.21(±0.03) + 0.34(±0.01)log(M em) /M*), r = 0.928). The first group comprises PMS 
stars with active accretion from the envelope. The third group comprises relatively highly 
extincted objects regardless of envelope mass. These are BMS stars with nearly edge-on 
disks (see Table [TJ. 

It is natural to expect a correlation between envelope mass and disk mass, but the 
results show no such correlation between these quantities. There is a strong concentration 
at (M env /M+) « io- LO±LO and (M disk /M+) « l(r 2 - 5±1 - 5 and there is a scattered population 
with (Mg^/M*) < 10~ 4 . The stars in the concentrated group are mostly stars in Spokes 
cluster or Cone(C). Their envelope mass of the scatter population may be an artifact of SED 
fitting as mentioned before. 

The envelope accretion rate is a strong function of envelope mass (log(M env /M^) = 
-4.21(±0.04) + 0.90(±0.02) • log(M env /M*), r = 0.944). Many active stars in the Spokes 
cluster or the Cone(C) region show an envelope accretion rate of io~ 5 - 8±0 - 5 M yr^ 1 . The 
tim e evolution of th e envel ope accretion rate obtained here is very similar to that in Figure 



4 of iRobitaille et al.l (120061 ). The outer radius of the envelopes does not show any correlation 



with age or mass. And there are only a handful of stars with R env ou t > 1000 AU. 



2.3. Median Properties 

We present the median values from SED fitter in Table [TJ Although the membership 
selection is far from complete in the Field region, we list the data for completeness. As the 
scatter of a given parameter is very large due partly to the age spread in a given SFR and 
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partly to the difference in environments, the values in the table are (logarithmic value of) 
median, 10 and 90 percentiles except the first two rows which give the fraction of objects in 
a given evolution stage. As shown in the 5th row of Table [T] the median age changes from 
0.35 Myr to 2.02 Myr. But if we define the age spread of a given SF R as the age diffe rence 
between the 10 and 90 percentiles in the cumulative age distribution ( iPark et al.l 120001 ). the 
age spread of SFRs in NGC 2264 is between 2.5 Myr and 3.4 Myr, which is larger than 
the difference in median ages among SFRs. This means that most properties of YSOs in a 
given SFR in parameter space overlap with those in other SFRs. Although we have such a 
limitation, we adopt the median age as the representative age of a given SFR. In this section 
we examine the variation of parameters among SFRs against their median ages. This picture 
can give an approximate view of several young open clusters with different ages. 



2.3.1. General Properties 



The first two rows of Table [J list t h e frac tion of objects at a given evolutionary stage. 
The definitions used in iRobitaille et al.l ( 120061 ) for the evolutionary stages are adopted and 
used here. The f raction of st a rs in Stages 0/1 or 2 is much higher than that of stars in Class 
I or II stages in ISung et al.l (120091 ). Such a difference is due mostly to the fact that only 
members or probable members (i.e. embedded objects in the three active SFRs) are taken 
into account. The difference in the definiti on of "class " and "stage" of YSOs also contributes 
to the difference. As already mentioned in lSung et al.l (120091 ) . the evolutionary stage of SFRs 
in NGC 2264 differs from one region to another, so we can expect difference in the age of 
each SFR. The age in Tabled] is the median, 10 and 90 percentile age of stars whose mass 
is between 0.2 - 1.0 M be cause most PMS evolutionary models give an abnormall y large 
age for massive PMS stars JSung et al.lll997l . I2OO0I . I2004T ) (see also iHartmann] (120031 ) I The 
Spokes cluster is the youngest {r me( nan = 0.35 Myr), S Mon, Cone(H), and Halo are slightly 
older (T me dian = 1-6 - 2.0 Myr), and Cone(C) is intermediate between them {r me dian = 0.76 
Myr). More discussions on the age of stars and star formation history in NGC 2264 are 
made in §3. 

To see whether there are any systematic variation in physical parameters among SFRs, 
we plot several parameters against the median age of a given SFR in Figure |5j The median 
value of interstellar extinction against median age is shown in Figure 0(1). Although we 
cannot find a noticeable trend for individual stars in a given SFR, the median relation shows 
a clear evolutionary trend between extinction Ay and age. Despite the small samples, young 
stars in the Field region or BMS stars deviate largely from the trend. Field members are far 
away from the active SFRs in NGC 2264 and so are more likely to be in relatively transparent 
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regions. The median Ay of BMS stars are slightly larger than the others. It is reasonable 
that BMS stars have a larger value of Ay dream as they are thought to have nearly edge-on 
disks. But the larger Ay of BMS stars may be just a result of small sample statistics or 
result from difficulties in fitting the SED of BMS stars. 

The median value of the total luminosity shows a strong correlation with age. As there 
is no correlation between age and mass, the correlation between age and total luminosity 
reflects the average evolutionary stage of a given SFR. A similar figure for individual stars 
in a given SFR shows a similar trend with a large scatter. Such a scatter is due to the 
diversity of stellar masses as well as to a difference in evolutionary stage. BMS stars show 
higher luminosity at a given age. That is probably due to the higher dust emission from 
disks around BMS stars as the median value of the disk masses of BMS stars is higher than 
for other stars. 



2.3.2. Disks 

Several of the next figures show the behavior of parameters relating to disk properties. 
The median value of the inclination angle of disks is the same for all SFRs and the distribution 
of inclination angle of individual stars for a given SFR does not show any preferential angle, 
i.e. a random distribution of inclination angles. On the contrary, the median value of the 
inclination angle of BMS stars is very high (the highest inclination angle from the SED 
fitter). This result is partly due to the difference in the best model selection criterion. If we 
use the same criterion for the best model for BMS stars, the median value of the inclination 
angle is 69.5° 4^5. A Kolmogorov-Smirnov (K-S) test between the cumulative distribution of 
inclination angles of PMS stars in the S Mon region and that of BMS stars gives 6.5 xl0~ 2 
% probability of a subset drawn from the same underlying distribution. This value is far 
lower than the probability between the S Mon region and the Spokes cluster (89.7 %) or the 
S Mon and Cone(C) regions (29.7 %). This result also implies that BMS stars are PMS stars 
with nearly edge-on disks. 

The mass in the disks as well as the mass accretion rate are strong functions of age, 
i.e. PMS stars in younger SFRs have larger disk masses and higher disk accretion rates as 
expected. The median value of the disk inner or outer radii is, on the other hand, smaller for 
younger SFRs and larger for older SFRs. The Rdisk,out of individual stars in an SFR show 
a similar trend as shown in Figure |3l but there is about two orders of magnitude difference 
between minimum and maximum value of Rdisk,out at a given age. The inner radius (Rdisk,in 
in unit of R su b or AU) of individual stars in an SFR does not show any trend as mentioned 
in §2.2.2. Other parameters, such as scale height factor Zf actor , flaring factor 0, and disk 
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accretion power a acc do not show any appreciable systematic difference or variation among 
SFRs. But the disk scale height at 100AU decreases as age increases. 

We also show the relation between disk mass and other parameters in Figure [5] (21) - 
(25). The disk accretion rate is a strong function of disk mass. The relation between the 
median value of the disk masses and the median value of the disk accretion rate is 

\ogM™% = -4.76(±0.35) + 1.28(±0.09) log , (r = 0.990) 

Disk mass and inner or outer radii are anti-correlated in the sense of median values. The 
intercept and exponent powers differ from the same relations for individual stars. The median 
value of the total luminosity is proportional to the median value of disk mass even for BMS 
stars, as seen in Figure [5] (24). But due to the variety of disk masses and evolutionary stages 
there is no correlation between disk mass and total luminosity for individual stars. The 
median value of the disk scale height at 100 AU is well correlated with the median value of 
the disk mass (Figure [5] (25)). This result means that although there is nearly no change 
in scale height factor Zf act0 r, the disk scale height at 100 AU is larger for massive disks, in 
other words it is larger at earlier stages. 

2.3.3. Envelopes 

The next figures (Figure [5] (13) - (18), (26) - (30)) show the time evolution of parameters 
relating to the envelopes of PMS stars in NGC 2264. It should be noted that in most cases 
the envelope could be observed for those stars with age < 1 Myr, therefore the median age 
of a given SFR cannot represent the age of individual stars with envelopes. 

The cavity opening angle increases as the median age increases. The decrease in median 
value of envelope mass with age is more dramatic being on average, a four orders of magnitude 
decrease in envelope mass from the youngest SFR Spokes cluster to the ~ 2 Myr SFRs S 
Mon, Cone(H) and Halo regions. The difference in envelope accretion rate among SFRs is 
small - only an approximate 5 times difference between the Spokes cluster and the S Mon or 
Halo regions. The envelope accretion rate with envelope mass of individual stars changes in 
a similar manner. The exponent of the power in the median relation (Figure [5)^26)) is 0.14, 
but that of individual stars is 0.85 

The median value of the envelope outer radii is well correlated with the median age of 
the SFR as shown in Figure [5] (14). This may be due to the infall of material to the disk. For 
individual stars there is almost no correlation between the envelope outer radius and stellar 
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age as mentioned in §2.2.3. The median value of the inner radius of the envelope increases 
as the median age increases as shown in Figure [5] (15). But any correlation between them 
is not so clear, because SED fitter gives the smaller, default value of the inner radius of the 
envelope (R en v,in = 1 Rsub) for many cases. 

The median value of the density in the cavity decreases with median age of the SFRs. 
But the median value of the ambient density p a mbient is nearly constant regardless of the age 
of the SFRs. BMS stars in Figure E^19) show a larger p am bicnt- As there is no reason to 
expect higher p am bient for BMS stars, this may be the result of difficulties of SED fitting for 
BMS stars. The circumstellar extinction (Ay from the outside edge of the YSO to the stellar 
surface along the line of sight) is a decreasing function of age as is the interstellar extinction 
in Figure 0T). The median value of the circumstellar reddening A V;Circum decreases as the 
median age increases with about the -4th power of the age. The Ay tCircum of BMS stars is 
very high due to extinction in the circumstellar disk, and does not appear in Figure [5](20) . 

Although there are strong correlations between the median value of envelope mass and 
the median value of disk mass or total luminosity as shown in Figures (29) and (30), there 
is no correlation between the same quantities for individual stars. 



3. STAR FORMATION HISTORY 



The star formation history in an SFR can give valuable information on star formation 
processes. From the age difference be tween massive stars near the MS turn-off and low- 



Herbia (I. 



962) proposed a sequential star formation 



mass PMS stars near the MS turn- on, 

scenario in an open cluster. But ISung et all ( 1l997l ) claimed that the mass-age relation 



of PMS stars in young open clusters may be an artifact of PMS evolution models, and 
suggested that such a discrepancy in age could be solved by r educing the PMS l ifetim e 
(Helmhoftz -Kelvin contraction time scale). Later ISung et all (120001) : iPark et al. (hoooh : 
Sung et al.1 (120041 ) supplied additional supporting results. However iPalla &: Stahler (2000) 
propo sed prolonged star formation which is nearly the same scenario as suggested by iHerbig 
fjl962h . 



Fiiresz. G. et al.1 ( 120061 ) identified three groups in NGC 2264 accord ing to their radia l 



velocity from high resolution multi-object spectroscopic observation. Later ISung et al.1 (120081 ) 
confir med the existence of t wo of them and found one more SFR, Halo, surrounding these 
SFRs. iTeixeira et al.1 (120061 ) first ident ified the Sp o kes cl uster as being dominated by Class 
I sources. Using Spitzer mid-IR data ISung et al.1 (120091 ) also identified two SFRs (Spokes 
cluster and Cone(C)) dominated by Class I objects and one SFR (S Mon) dominated by 
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Class II objects. Throughout this paper we use the name of the SFRs used by ISung et al. 
( 120091 ). It is very interesting to now compare the age distribution among SFRs. It should 
also give some important information relating to the star formation processes within NGC 
2264. 



3.1. Age Distribution 

The age of PMS stars is determined in two ways. One is the age from the SED fitter, and 
the other is the age from the CMDs (see §2.1.2). We present the cumulative age distribution 
of low-mass PMS stars in a given SFR in Figure [6J To avoid unrealist ic ages from PMS 



evolu tion models, we limit the PMS stars with m = 1.0 - 0.2 M (see ISung et al.l ( 11997I . 



2004m . 



We also present the median age, 10 and 90 percentile in the cumulative age distribu- 
tion, and the results from K-S tests in Table [2J If we assume that the median age is the 
representative age of an SFR and the difference between the 10 and 90 percentile in the age 
distribution is the age spread of the SFR, the median age from SED fitter is about 1-2 Myr 
younger than that from CMDs for a given SFR. In addition, the age spread is about 3 Myr 
from SED fitter, while that from CMDs is about 5-6 Myr. The smaller value of age spread 
is due to the fact that the stars used in SED fitter are relatively younger and easier to detect 
with Spitzer observations. On the other hand, the somewhat larger age spread from CMDs 
may be due to the effect of reddening. In addition, the reddening law in young SFRs may 
differ from that in the general interstellar medium. Therefore, the age spread, or in other 
words the cluster formation time scale, may be about 5 Myr. 

The left panel of Figure [6] shows the cumulative age distribution of PMS stars from 
SED fitter. Although the median age of stars in Field is much younger than that in Halo, 
this seems to be due to the small number of member stars in Field. As the age distribution 
of stars in Field is very similar to that of Halo and the number of known members with 
Ha emission or X-ray emission in Field is very small, we merge the stars in Field and Halo. 
The age distribution of stars in a given SFR is far different from each other. The stars in 
the Spokes cluster are the youngest among the SFRs in NGC 2264. The stars in S Mon 
and Cone(H) have a similar distribution, while the age of stars in Cone(C) is intermediate 
between the Spokes cluster and S Mon. The age of stars in Halo is slightly older than that in 
S Mon, but the difference is marginal. To check whether the age distribution among SFRs 
in NGC 2264 arise from the same distribution or not, we have performed a K-S test. The 
results are listed in Table [2j The probability of the subsets being drawn from a common 
population is very low between S Mon and the Spokes cluster or S Mon and Cone(C), but 
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is very high between S Mon and Cone(H). The probability between S Mon and Halo (and 
Field) is marginal. 

The cumulative distribution of ages from the CMDs is presented in the right panel of 
Figure [61 The pattern of age distributions is far different to that from SED fitter. The age 
distributions of all active SFRs (S Mon, Spokes cluster, Cone(C), and Cone(H)) in NGC 
2264 are very similar. As many young PMS stars in the Spokes cluster or Cone(C) are 
deeply embedded, they are not in the PMS locus of NGC 2264. In addition, the PMS stars 
in the PMS locus are situated in the outside surface region of NGC 2264, and therefore 
the contamination of Halo or Field members is unavoidable. This is the reason why the 
age distribution of these active SFRs is similar to each other. On the other hand the 
difference in age distribution between stars in S Mon and Halo (including Field) is rather 
clear. The probability from a common population is very low (P = 4 x 1(T 4 ). Somewhat 
larger probabilities relative to the age distribution of stars in the Spokes cluster are obtained 
from the K-S test. The results are due to the small number of member stars in the PMS 
locus of NGC 2264 because most of the members of the Spokes cluster are an embedded 
population. In order to check the reliability of this comparison, although there is no way to 
estimate the age of non-member stars, we estimated the age of non-member stars in the PMS 
locus of NGC 2264 by assuming they were members of NGC 2264. The age distribution of 
non-member stars is far different from that of Halo or other member stars in active SFRs. 
The very low probabilities from a K-S test between the age distribution of non-member stars 
in Field and other SFRs indicate that this kind of test is valid for the study of star formation 
history. 



3.2. Triggered Star Formation in NGC 2264 



The age distributions obtained above give a clear picture of sequential star formation 
in NGC 2264. The stars in Halo (and Field) region are formed first, followed by stars in 
S Mon and Cone(H). The stars in Cone(C) formed later and finally star formation took 
place in the Sp okes cluster. We can find several Clas s I objects in Cone (C), and many 
Class I objects (jSung et al.l 120091 ) and Class objects (ITeixeira et al.l 120071 ) in the Spokes 
cluster. In other words, star formation activity in NGC 2264 started at the surface of the 
giant Monoceroti s molecular cloud and propagated inside the cloud. As already noted in 
Sung et al.l (120081 ) there are several small SFRs spread over Field region. The median age 
of Field is very similar to that of Halo. It is very difficult to imagine that spontaneous star 
formation started first in a low density region (Halo or Field) of NGC 2264. If the star 
formation history derived in §3.1 is true, we should accept that star formation in NGC 2264 
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was triggered by an external source, such as a supernova (SN) explosion that occurred about 
6-7 Myr ago. 



An extensive review of triggered star formation has been published by lEfmegreenl (119981 ). 
He classified trigger ed star formation into three types according to the scale of star formation. 
Tauber et all ( 119931 ) conducted a high resolution CO observation in the bright nebular region 
of NGC 2264 and found three bright clumps. This is an example of small scale triggered star 
formation in NGC 2264. The age difference between S Mon and Cone(C) is 0.4 - 0.8 Myr, 
and the projected distance between S Monocerotis and HD 47887 is about 26' (= 5.7 pc at 
760 pc). If star formation in Cone(C) were triggered by the influence of S Monocerotis, the 
minimum propagation velocity of star formation activity would have been 7-14 km s _1 , 
which is far faster than sound speed in a general molecular cloud (~ 1 km s^ 1 ). Therefore 
intermediate scale triggering - "collect and collapse" - may not have taken place in NGC 
2264. 

The first mentio n of large sca le triggered star formation happening in and around NGC 
2264 is the pap er by iHeilesI (119791 ) who found an HI shell in the Monocerotis region. Later 
Schwartzl (119871 ) found a half circle stellar loop in Mon OBI region using IRAS images. He 
listed several young objects along the semicircle - IC 446, IC 2169, NGC 2264 IRS1 and 
IRS2, and NGC 2261 (= R Mon). The dia meter of the half circle was about 2.5 degree 
(about 33 pc at 760 pc). lOliver et al. Jl996h studied the large scale structure in Mon OBI 
from an extensive CO survey. They identified many molecular clouds and listed candidate 
optical objects possibly associ ated with them. T he center of the shell is (I, b) — ( ~ 202°, 
~ 2°) from Figure 1 and 2b of lOliver et al.l (119961 ) . 



3.3. Objects of Common Origin with NGC 2264 

We showed in the previous section that the star formation in NGC 2264 was triggered 
exte rnally, and the t r iggering source ma y have been an SN that exploded about 6-7 Myr 
ago. ISchwartzl (119871 ) : lOliver et all (119961 ) listed many possible optical counterparts. We have 



checked these objects and a few should be excluded. All data were obtained from the Simbad 
database |f| and the open cluster database WEBD A.B There ar e no op tical counterparts for 
approaching clouds 4, 12, and 19 in Figure 2a of Oliver et all ( 19961 ). Although the large 
extent in velocity of cloud 4 in Figure 1 of lOliver et al.l (119961 ) implies the possibility of it 



http : / /simbad . u-strasbg . f r 



http : //www . univie . ac . at/webda 
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being an ejected cloud from the imputed SN explosion, the cloud is located far from the 
center of the CO hole (I ~ 202°, b w 2°). Cloud 4 may therefore not have the common 
origin. Cloud 19 may also not be related to the SN explosion for the same reason. Cloud 9 
shows a large velocity gradient and is very close to cloud 10. There are two reflection nebulae 
(IC 2167 = IC 446; GN 06.28.4) in cloud 9, and several young variable stars in and around 
these reflection nebulae. These young stars are 2MASS J06310611+1027340 (Sp: B2.5V) in 
GN 06.28.4, VY Mon (Sp: B8, Ha emission star) in IC 2167, V540 Mon (= HBC 521, Sp: 
M6, Ha emission star) and V687 Mon (= LkHa 274, Sp: K4IVe) around IC 2167. 

Cloud 10 is very close to the CO hole and is an approaching cloud. Within it are an 
HII region IC 2169 and 4 early type stars - HIP 31038 (= vdB 76, Sp: B7IIIp), V727 Mon 
(Sp: B8V, eclipsing binary), HD 258749 (= vdB 77, Sp: B5III), and HD 258853 (= vdB 
78, Sp: B3). There are two reflection nebulae in cloud 11 - NGC 2245 and NGC 2247. An 
Ha emission star V699 Mon (= LkHa 215, Sp: B6) is at the center of NGC 2245, and a 
B6pec star V700 Mon is in NGC 2247. Although cloud 11 is a very massive cloud, only 
two reflection nebulae indicate current star formation activity. This may be due to its small 
cross section facing the CO hole. 

Cloud 13 is related to IRAS 06337+1051, many radio sources and a thin nebula around 
the mid-IR source. Several young or massive star candidates can also be found. Among them 
Kiso Ha 39 may be related to the cloud, but the rel ation of GSC 00 7 37-01170 (Sp : 05:) 
and Tycho 737-898- 1@ (Sp: B) to cloud 13 is uncertain. ISchwartzl Jl987h : bliver et all Jl996h 
proposed a possible relation between cloud 13 and the open cluster NGC 2259. However, 
NGC 2259 is an intermediate-age ope n cluster that canno t be related to the cloud. The most 
massive molecular cloud in Table 1 of lOliver et all (119961 ) is the giant molecular cloud which 
is associated with NGC 2264. A small cloud - cloud 17 - harbors a reflection nebula NGC 
2261 and the young star R Mon (Sp: B0). 



Oliver et all (119961 ) proposed a possible relationship between the open cluster NGC 2254 



and clouds 18 through 20. NGC 2254, as well as Basel 7 and Basel 8, have been observed only 
in the photographic RGU system and so their nature is somewhat uncertain. They appear 
to be intermediate-age op en clusters from th eir CMDs. Another open cluster NGC 2251 
is marked in Figure 2b of lOliver et all (I1996I ). The cluster's distance in WEBDA is more 
distant than NGC 2264, and although the data were derived mostly from UBV photographic 
photometry the cluster seems to be a Pleiades-age open cluster. Therefore the relationship 



5 V = 10.18, B - V = 1.22, U-B = 0.42 bv lDrillind (|1975[ ). The B — V color is too red for a B type star. 
If the star were a low-mass PMS star and its blue color caused by UV excess, the star would be somewhat 
brighter than the low-mass members of NGC 2264. 
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between the clusters and the clouds is quite uncertain. New homogeneous data are required 
for a comprehensive view of the star formation history in the Mon OBI association. 



THE INITIAL MASS FUNCTION 



In this section we derive the IMF of NGC 2264 in two ways - using masses derived 
from SED fitter and masses derive d from fitting CM Ds. The membership selection is nearly 
complete do wn to about 0.25 M (ISung et al.ll2004l ) for two Chandra field s - one around S 



Monocerotis (iRamfrez et al.l 120041 ) and the other around the Cone nebula (jFlaccomio et al. 



20061 ). These two fields fully cover the four active SFRs (S Mon, Cone(C), Cone(H), and 
Spokes cluster) in NGC 2264. We designated these four regions as the reference region for 
the combined IMF. The IMF is defined as the number of stars per unit projected area per 
unit logarithmic mass interval, i.e. £(logm) = dN/ A log m /area , where area is in kpc 2 . We 
adopted the distance of NGC 2264 as 760 pc (jSung et al.lll997l ). We used the logarithmic 
mass interval Alogm = 0.1 for the reference region, and Alogm = 0.2 for the IMF of an 
individual SFR. We also calculated the IMF for each SFR in the same bin size, but shifted 
by 0.1 in logm to reduce the binning effect. For massive stars (logm > 0.5) there were 
many cases with no star in a given mass bin. In those cases, we enlarged the mass interval 
so that there was at least one star in each mass bin. The upper mass limit was assumed to 
be 100 M - the most massive star ever measured was 116 ± 31 M (ISchnurr et al.ll2008l ) - 
and we calculated the mean value of log m in the mass bins. One thing to be noted is that 



although some st ars were known to be in a binary system, such as S Mon (IGies et al.l Il996 



Sung et al.lll997l ). we assumed all stars to be single as we have only limited information on 
multiplicity (see §6.1 for the effect of binarity). 



4.1. The IMF from SED Fitter 

We determined the IMF for each SFR, but show the IMF of the reference region in 
the left panel of Figure [TJ because a) membership selection is nearly complete; b) we could 
determine the slope of the IMF in the massive regime (M > 3 M ) reliably as only a few 
massive stars were in a given SFR. The right panel of Figure [7J shows the relative IMF of 
each region relative to the reference IMF. The abrupt decrease in the IMF at logm < -1.0 is 
due to the limited exposure time of the Spitzer observations of NGC 2264. There are some 
stars with logm < -1.0 although the lowest mass from SED fitter is 0.1 M ; the mass of 
these stars is derived from CMDs only. 
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The IMFs by iKroupal (120021 ) and by IMoraux et al.l (120031 ) are superimposed for com- 
parison. Their IMF were adjusted roughly at logm = -0.3 - -0.4. There is a distinct peak in 
the combined IMF of NGC 2264 at about 2 M and another less distinct peak at about 0.4 
M as already mentioned i n §2.1 . 2. We could not find similar peaks in the IMF of field stars 



young open clusters (jSung et al 



in the Solar neighbor hood dScalo Il986l ; iLee fc Sunall995l; IChabrierll200ll ; lKroupall2002l ) or in 



2004t iMuench et al.ll2002l ). In many cases there was a local 



maximum at about 1 M in the IMF of field stars or open clusters, therefore the prominent 
peak at M 2M in Figure UJ may be an artifact induced by the \ 2 minimizing routine in 
SED fitter. Such a peak can be seen the YSO ma ss function of M17 SWex (see §6.2 for the 
discussion on this issue) (jPovich &: Whitney! |2010[ ). Normally the IMF can be characterized 
by its slope T(= d log (,/d logm). The slope fo r massive stars (lo gm > 0.5) is -1.79 ± 0.14, 
which is very similar to the slope obtained by iPark et al.l (120001 ) for NGC 2264 or by IScalo 
( 119861 ) for the Solar neighborhood. 

Figure [7J (right panel) shows the difference in the IMF for each SFR relative to the 
combined IMF in the left panel. The IMF itself, by definition, represents the surface density. 
The surface density of stars in the Spokes cluster or the Cone(C) region is higher (Alog£ = 
+0.43 ± 0.12 in Spokes cluster and Alog£ = +0.24 ± 0.15 in Cone(C) in the mass interval 
of logm = +0.5 - -0.5) than the average value of four active SFRs. And Cone(H) is slightly 
lower (Alog£ = —0.31 ± 0.18 in the same mass interval). The IMF of the Spokes cluster 
shows a slight enhancement at about 2.5 M above the mean (about 2 a level). The surface 
density of the Halo region is very low (A log£ = —0.76 ± 0.17) and lacks massive stars. The 
low surface density may be due partly to the incompleteness of the membership selection 
in the Halo region. The IMF of the Field region is far lower. This is due partly to the 
incompleteness of membership selection but largely to the huge area observed relative to the 
small SFRs spreading over the Field region. 



4.2. The IMF from Color-Magnitude Diagrams 



The IMF of open clusters is derived from their CMDs. But for young open clusters 
where the distribution of reddening material is very inhomogeneous, there are many limita- 
tions, especially for low-mass PMS stars, because there is practically no ways to determine 
the reddening of individual stars from photometric data alone. In addition, there is the un- 
certai nty in PMS evolu tion models. We derived the IMF of stars in the PMS locus of NGC 
2264 ( ISung et al.ll2008l ) by applying the mean reddening. In this case we could not take the 
heavily embedded stars and BMS stars into account. We derived the IM F of NGC 2264 fo r 
two values of reddening - the mean reddening (E(5 — V) — 0.07 mag - ISung et al.l ( 119971 ); 
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Park et all ( 120001 ) ) determined from early type stars in NGC 22 6 4, an d the median Ay from 
SED fitter for each region. As already checked in ISung et al.l (120041 ) a small difference in 
reddening does not greatly affect the shape of the IMF, i.e. the more important factor is 
the inclusion or exclusion of embedded stars and BMS stars. As PMS stars in the PMS 
locus are mostly less reddened, a small value of E(B — V) is more appropriate. The IMF 
determined from the stars in the PMS locus is shown in Figure The IMF presented here 
is the combined IMF of the reference region. The only massive star (M > 10 M ) missed in 
the derived IMF is NGC 2264 IRS1, an embedded flat SED object with 13.5 M G from SED 
fitter. The shaded area represents the mass range where membership selection is largely 
incomplete as mentioned before. 



The IMF of NGC 2264 is very similar to that of the Pleiades (IMoraux et al.ll2003[ ) down 
to logm ~ — 1.2 where the membership selection is incomplete even in the reference region 
of NGC 2264. But there is some scatter for massive stars (m > 2 M Q ). This is related to the 
fluctuations due to the small number of massive stars in the cluster. The IMF of NGC 2264 



shows a slight deficiency in massive stars relative to the IMF of iKroupal (120021 ) . We could 
not find a similar peak in the IMF at about 2 M . The sl ope of th e IMF is r = —1.6 5 ±0.1 4 
for logm > 0.5. This slope is much s teepe r than that of ISalpeterl (119551 ) or IKroupal (120021 ) 



but is very similar to that of IScalol (119861 ) ; lHambly et al. Jl999h . In this case the surface 



density in the Spokes cluster is very similar to that of S Mon because many members of the 
Spokes cluster are heavily embedded. The surface density of Cone(C) is slightly elevated, 
and that of Cone(H) is slightly lower. 



YOUNG BROWN DWARF CANDIDATES 



One of main purposes of this study was to select the YBDCs in NGC 2264. NGC 2264 is 
one of the nearest young open clusters and the foreground extinction is practically zero. The 
YBDCs in NGC 2264 are therefore still very b right and their photo metric characteristics are 
very similar to those of young low-mass stars. [Kendall et al.l (120051 ) tried to find the YBDCs 
in NGC 2264 using deep Iz photometry with CFH12K. Although they tried to find YBDCs 
in various color-color diagrams and CMDs and through the help of recent model atmospheres, 
they were unable to find YBDCs confidently due to a lack of reliable membership selection 
criteria. They presented eight probable YBDCs and two additional candidates. 

We selected YBDCs with the following criteria. 



(1) the mass from CMDs is smaller than 0.075 M Q . As ISiess et al.l (120001 ) published 
PMS evolution models only for m PMS > 0.1 M Q , we have to use the PMS evolution mod- 
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els by IBaraffe et al.l ( 19981) for th e selec tion of YBDCs. But as it is well known that the 
theoretical colors of IBaraffe et al.l ( 119981 ) are not well matched for low-mass stars, we have 
transformed their model parameters {T e ff and Lr„ ; ) to V — I or R — I using the empirical 
color-tem perature and color-B C relations by iBesselll (ll99lf ). Still the theoret i cal ev olutionary 
tracks by IBaraffe et al.l ( 119981 ) do not match well with those of ISiess et al.l (120001 ) . but it is 
no w possible t o extr apolate the tracks to the red color regime by following the 0.1 M track 
of ISiess et all (boooh . 



(2) Stars with one or more membership criteria such as Ha emission, X-ray emission, 
or IR excess (Class I, II, or III objects). 

(3) The mass from SED fitter if the star is detected with the Spitzer should be 0.1 M 
which is the lowest mass from SED fitter. 

Using these criteria, we have selected the 79 YBDCs presented in Table [3j The CMDs 
of these YBCDs are shown in Figure [9j The PMS evolutionary tracks are superimposed in 
the figure. IR excess YBDCs are relatively bright due to the limited exposure time used in 
the Spitzer observations. In addition, X-ray emission cannot identify YBDCs fainter than 
/ « 19 and while Ha photometr y reaches down to fainter limits, the selection probability is 
somewhat low (about 30% - see ISung et al.l ( 12004 ) ). 



Three of the YBDCs in Table [3] (C33215 = 2MASS J06405674+0938101, C35463 = 
J06410604+0949232, and C20186 = J06401789+0941546) are in common with the YBDCs 



listed bv Kendall et al.l (120051 ). Another object (2MASS J06402759+0945464) could be a 
YBCD because the star is within the PMS locus of both CMDs and R— Ha is very similar to 



that o f weak line T Tauri stars with X-ray emission. The other 6 YBDCs in Kendall et al. 



(120051 ) may not be members of NGC 2264 - Three objects (2MASS J06395722+0941011, 
J06400642+0 944197, J0641128 1+0945529) are fainter than the lower boundary of PMS locus 
of NGC 2264 rtSung et al.ll2008h . Two (J06401053+0939557, J06413132+0935120) are within 
the PMS locus in the (J, V — I) diagram, but are below the lower boundary of the PMS 
locus in the (I, R — I) diagram. In addition the R— Ha colors of these two objects are very 
similar to those of reddened background stars. And 2MASS J06404873+0939017 is detected 
only in /. 
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6. Discussion 

6.1. The Effect of Binarity on the IMF 

The fraction of binary or multiple systems and their ch aracteristics are impo rtant con- 
straints for the studies of stellar formation and evolution ( IBurgasser et al.l 120071 ). Several 
studies on the frequency of multiple systems and their properties were made especially for 
very low- mass stars (see ISterzik fc Durisenl (120031 ) and references therein), but their initial 
values are still not well known. The frequency and mass ratio in multiple systems may affect 
the actual shape of the IMF. 

In general the IMF derived from young open clusters is just the IMF of single stars and 
primaries of multiple systems. It is very important to check the effect of multiple systems 
on the IMF. We constructed a model cluster of 1,000 primaries plus single stars with a 
Monte Carlo method. The IMF of NGC 2264 in Figure E is adopted as an input IMF. The 
frequency of binary as a function of prima r y mas s was adopted from the results of two- 
step dynamical model of ISterzik fc Durisenl (120031 ) which reproduce well the observations. 
In addition the secondary t o primary mass ratio as a function of spectral type in Figure 4 
of ISterzik fc Durisenl (120031 ) was converted to a function of primary mass and used in the 
simulation. We assume that the binary fraction of massive stars (m > 10 M ) was 90%. As 
there is nearly no information on the fraction and mass ratio in multiple systems of triple or 
more, only single stars and binary systems were taken into account. 

A total of 355 binary systems and 645 single stars were generated (see Figure [TUT) . The 
massive part of the mass function is dominated by the primaries in binary systems, while 
the low-mass part is due to single stars. The mass function of secondaries is very similar 
to the IMF of all stars. The shape and slope at the massive part of the IMF of all stars is 
very similar to those of the input IMF. In addition, the peak of the IMF occurs at logm « 
-0.6, which is the same value as the input IMF. Such similarities are due to the facts that 
the binary fraction of very low-mass stars and BDs is nearly zero and that of massive stars 
is very high. In addition the mass ratio of very low-mass stars and BDs is relatively high, 
but that of massive stars is preferentially a smaller value. 

If we change the fraction of binary systems of very low-mass stars and BDs from 0% to 
10 %, the fraction of binary systems increases from 35.5 % to 37.1 % (371 binary systems 
were generated in this case). Although the mass of the IMF peak does not change, the slope 
of the IMF of the massive part becomes steeper and the IMF of very low-mass and BD 
regime stars increases. 



We can conclude that if the binary fraction and mass ratio distribution are very similar 
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to those adopted here, the shape and peak of the IMF do not change very much. But the 
fraction of multiple systems and the mass ratio distribution can affect strongly the IMF of 
all stars. And therefore it is very important to investigate the frequency and mass ratio 
distribution of multiple systems for several open clusters with different ages. 



6.2. YSO Mass Function of M17 SWex 



Recently iPovich &: Whitney! (120101 ) studied the YSO mass function (YMF) of the em- 



bedded massive star forming region Ml 7 SWex. They used SED fitter to e stimate the mas s 



of YSOs in M17 SWex and found the YMF to be much steeper than that of ISalpeterl ( 119551 ). 
Their YMF showed a peak at m ~ 3M . The location of the YMF peak is very similar to 
the peak in the IMF in Figure [7J The prominent peak in the YMF of M17 SWex may be 
due to the combination of two effects - the completeness of the PMS membership selection 
in M17 SWex and the effect mentioned in §4.1. 

If we attribute this peak to an artifact from SED fitter, then the slope of the YMF of 
Ml 7 SWex may be somewhat shallower. In addition, if we take account of the fact that 
there are no known very massive young stars (e.g. M > 3OM ) with disks, the most massive 
stars in the region may have been missed and the slope of the YMF of Ml 7 SWex could be 
normal. 



7. SUMMARY AND CONCLUSION 

We have studied the young open cluster NGC 2264 using all available photometric data 
fro m optical to mid - IR. W e performed SED fitting using the \ 2 minimizing SED fitting tool 



by iRobitaille et al.l ( 120071 ) for more than 1,000 stars detected in all Spitzer IRAC channels 
with high quality (ejRAc,totai < 0.25 mag) mid-IR data. We estimated the mass and age 
of member stars in NGC 2264 and constructed the star formation history and IMF of the 
cluster. The results obtained are summarized as follows. 

(1) The parameters relating to the central stars from SED fitter and other methods are 
compared. The temperature from SED fitter is, in many cases, higher than that expected 
from the spectral type of the star. And the mass of central stars shows a conspicuous clump 
at logm ~ 0.2 - 0.3 and -0.2 - -0.4, which causes anomalous peaks in the IMF from SED 
fitter. The age of central stars is, on the other hand, smaller for sub-solar mass stars, but is 
larger for massive stars (m > 1 M ). 



(2) The parameters relating to disks or envelopes of individual stars in NGC 2264 were 
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examined and found to be very similar to those in iRobitaille et al.l (120061 ) . The properties 
of individual stars showed no significant differences among the SFRs in NGC 2264. As the 
maximum difference in median age among SFRs is smaller than the age spread in an SFR, 
the large scatter of parameters is due mainly to their age spread. 

(3) The median values of the physical and structural parameters showed an evident 
difference among SFRs in NGC 2264 that was strongly related to the median age of the 
SFRs. Such a difference among SFRs gives an approximated view of several young open 
clusters with different ages. 

(4) The cumulative age distribution of stars showed distinct differences among SFRs. 
The median age of SFRs implies that star formation in NGC 2264 occurred sequentially - star 
formation started at the surface of the Monocerotis giant molecular cloud and propagated 
into the cloud, i.e. the stars in the Halo and Field region are the oldest, the stars in S Mon 
and Cone(H) are intermediate while those in Cone(C) are younger, and those in the Spokes 
cluster are the youngest. The star formation history seen in NGC 2264 and the fact that the 
first star formation occurred in the low density region (Halo and Field) of NGC 2264 imply 
that star formation in NGC 2264 was triggered externally. 

(5) The IMF using the stellar mass from SED fitter shows a prominent peak at m « 
2M . As such a peak cannot be found in the IMF of other open clusters or field stars in the 
Solar neighborhood, the peak in the IMF from SED fitter may be an artifact of the SED 
fitting tool. The IMF from CMDs is well consistent with the IMF of the Pleiades. 

(6) The slope of the IMF of massive stars (logm > .5) is -1.7 ± 0.1. This value is 
somewhat st eeper t han t he standard Salpeter-Kroupa IMF (jKroupall2002l ). but is consistent 
with that of lScalol (jl986h . 



(7) We selected 79 young brown dwarf candidates from the PMS loci in the CMDs of 
NGC 2264 that showed either Ha emission, X-ray emission, or IR excess emission, and with 
implied masses smaller than 0.075 M . 
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Fig. 1. — Spectral type versus temperature relation. (Left panel) The effective temperature 
from SED fitting. (Right panel ) The ef f ective temperature is derived from the V — I or R — I 
versus temperature relation of iBesselll (1l995[ ). We applied the same red dening for all stars 
(E(B — V) — 0.07). Spectral types are taken from iRebull et al.l ( 120021 ). Dots, circles, and 
crosses represent, respectively, Class II objects, stars with normal photosphere, and stars 
without X-ra y and Ha emission. T he solid line denotes the spectral type versus temperature 
relation from ISchmidt-Kalerl ( 119821 ). 



31 




log m(M ) CMD age CMD (Myr) 



Fig. 2. — Comparison of stellar mass and age. (Left panel) Comparison of stellar mass from 
color-magnitude diagrams with that from SED fitter. The size of symbols is proportional 
to the reddening from SED fitter. Dots and open squares represent Class II objects and 
stars with normal photospheres, respectively. (Right panel) Comparison of age from color- 
magnitude diagrams with that from SED fitter. Circles and squares represent Class II objects 
and stars with normal photospheres, respectively. Filled and open symbols denote objects 
with M^sed < 1 M and M^sed > 1 M , respectively. 
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Fig. 3. — Disk parameters of individual stars from SED fitter. The upper two panels show 
the time evolution of several parameters relating to disks around PMS stars, while the lower 
panel displays the variation of several parameters against disk mass. To see whether there 
are any systematic differences among SFRs or evolutionary stages we use differen t colors and 



differ ent symbols. Dots and crosses represent, respectively, Class I or II objects (jSung et al. 



20091 ) and stellar photospheres. Blue, black, red, green, yellow, and cyan represents the PMS 



stars in S Mon, Spokes cluster, Cone(C), Cone(H), Halo, and Field, respectively. 
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Fig. 4. — Envelope parameters of individual stars from SED fitting. The upper panel shows 
the time evolution of several parameters relating to the envelopes. The lower panel shows 
the relation between circumstellar reddening, disk mass, or envelope accretion rate and the 
mass in the envelope. Symbols are the same as in Figure [H 
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Fig. 5. — Results from SED fitting. To see the evolutionary trend of YSO parameters, we 
show the median value of young stars in a given SFR. A blue dot, black square, red star, green 
asterisk, triangle, small open square and small cross represents the median value of young 
stars in S Mon, Spokes cluster, Cone(C), Cone(H), Halo, Field, and BMS stars, respectively. 
The actual scattering of individual stars in a given SFR is very large and therefore in most 
cases, it is very difficult to see the evolutionary sequence for the region. In addition, the 
relations for envelopes are only for small samples of stars with envelopes in a given SFR. See 
text for details. 
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Fig. 6. — Star formation history of NGC 2264. (Left panel) Age of PMS stars from SED 
fitter. (Right panel) Ages from CMDs. We assumed E(B — V) = 0.07. To avoid unrealistic 
ages from PMS evolution models, we limit the ages to PMS stars with m = 1 - 0.2 M Q . 




Fig. 7. — The IMF of NGC 2264 from Syitzer obs ervations. (Le ft panel) The IMF of stars in 



the S MON and CONE regions (large dots) of lSung et al.l (12008 ) wher e membership selection 

(120031 ) for the Pleiades (open 



Moraux et al. 



is nearly complete d own to 0.25 M^ . The IMF of * 

circle) and the IMF iKroupal (120021 ) (dashed line) are superimposed for comparison. (Right 
panel) Differences in the IMF of each region relative to the IMF in the left panel. 
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Fig. 8.— The IMF of PMS stars in the PMS locus of NGC 2264 (large dots). We assumed 
E(B - V) = 0.0 7. (Left panel) The IMF of stars in the S MON and CONE regions of 



Sung et al.l (120081 ) where membership selection is nearly complete down to 0.25 M Q . (Right 
panel) Differences in the IMF of each region relative to the IMF in the left panel. 




Fig. 9. — Color- Magnitude diagrams of you ng brown dwarf c andidates in NGC 2264. Small 
dots are all stars in the field of NGC 2264 (ISung et al.ll2008l ). The R — I color s of ve ry red 
stars were corrected using the prescription described in Appendix of ISung et al.l ( 120081 ). The 
thick dashed lines represent the locu s of PMS stars in NG C 2264. Thin solid an d dotted lines 
show the PMS evolution tracks from lsiess et al. ( 2000 ) and Baraffe et al. ( 1998), respectively . 



The empirical color-temperature and color-BC relations for very red stars from 
were used. 



Besselll (I199lh 
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Fig. 10. — The effect of binarity on the IMF. A Monte Carlo simulation was performed to 
construct a model cluster of 1,000 primary plus single stars. The input IMF is the IMF 
of NGC 2264 in Figure El The binary fractio n and the mass rat i o bet ween primary and 
secondary as a function of primary mass is from lSterzik fc Durisenl (120031 ) . Square, triangle, 
and star symbol represent, respectively, the IMF of primaries, secondaries, and single stars, 
while large dots stand for the IMF of all stars including the secondaries of binary systems. 
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(145) 


-0.09 


0.00 
-.15 


(100) 


-0.09 


0.00 
-.15 


(177) 


-0.11 


-.02 
-.15 


(33) 


-0.27 


1.02 
-.92 


(257) 


0.12 Z% 


(176) 


-0.07 


1.18 
-.67 


(145) 


-0.36 


0.96 
-.79 


(100) 


-0.40 


0.99 
-.86 


(177) 


-0.07 


1.70 
-.92 


(33) 



56.6 


Hi (145) 


56.6 


Hi (100) 


56.6 


m (177) 


56.6 


-3.32 


ZIU (145) 


-4.84 


l7 97 (100) 


-4.29 


:?:s4 (177) 


-3.88 


2.19 


?.|(145) 


2.40 


?:ii (100) 


2.37 


lit (177) 


2.11 


1.51 


0:32 (82) 


2.11 


IM (72) 


2.13 


Ut (126) 


2.09 


-0.58 


-? 8 44 (145) 


0.43 ^ 3 48 (100) 


0.00 


T« (177) 


-0.96 


-8.89 


~-% 8 L (145) 


-10.93 


:i 4 4 36 (100) 


-10.10 


Ii3 8 67 (177) 


-9.86 


0.80 


8:1 (145) 


0.85 


8:1 (100) 


0.82 


8:11 (177) 


0.78 


1.10 


lit (145) 


1.10 


lit (100) 


1.09 


lit (177) 


1.10 


-1.80 


^2.77 (145) 


-1.91 


zlil (100) 


-1.86 


"JJo (177) 


-1.79 


4.91 


2° 7 8 (145) 


4.43 


!: 3 76 (100) 


4.49 


III (177) 


5.06 



-21.; 



-21.2 
-22.0 



(145) 

-1.70 Zi\ 2 (145) 



-21.8 I220 (10°) 
-3.03 °f 09 (100) 



-21.8 



_ 2 07 0.76 
z - y ' -4.10 



IH (177) 
(177) 



(33) 



!|:8 4 (33) 

-7.43 (33) 
2.99 
1.66 

Hi (14) 

-? 8 53 (33) 
=I 3 5 30 (33) 

? : 1 (33) 

1:04 (33) 

-2.81 (33) 
8.84 (33) 



3.31 

46.3 
21.3 

3.08 3)01 

2.54 l-H 



30.6 



-5.00 



(14) 
(14) 
(5) 
(33) 



-irsi 
-8.20 

-6.87 ZIU (14) 
-20.7 zfH (16) 



-91 Q 21 - 2 
z± - a 22.0 

2 30 2 - 43 



(33) 
(33) 



1 



0.10 
0.85 
•41 Ht (19) 
1-46 HI 
-0.10 (40) 
0.21 L- 4 5 °4 (40) 



87.1 

-2.79 
2.36 
1.16 

-0.51 

-9.12 
0.78 
1.08 

-1.76 
4.23 

49.7 : 



Hi (40) 
'-All (40) 
HI (40) 

8:i (22) 

H 8 25 (40) 
="33 (40) 

8:11(40) 

lit (40) ' 

=i:il (4o)w 
5:il (40) 



(16) 



31.9 

3.38 3 -^ (16) 



0.86 
-5.36 
-6.90 



3.03 

lil (10) 
(40) 



-7.60 



-21.5 



-21.0 
22.0 

2.19 J 2 .78 



(16) 



(40) 
(40) 



a The data given are median, 10%, 90%, and number used in the calculation. 
b We choose the Xmin model as the best model for BMS stars. 

d Due to the uncertainty in the age of massive PMS stars, the mass limit used here is m = 0.2 - 1.0 Mq. 
d We neglect the data whose value is 1.00 R su b- 



Table 2. Ages and K-S Test Results 



method 


SFR 


N st ar a 


Tmedian (Myr) a 


no% (Myr) a 


T90% (Myr) a 






Probability from Kolmogorov- 


Smirnov Test 




reference 




O 1V_HJI1 




Spokes cluster 


Cone(C) 




Halo 




S Mon 


120 


1.63 


0.27 


3.27 
















Spokes cluster 


78 


0.35 


0.02 


2.55 


1.65x10" 


8 










SED 


Cone(C) 


82 


0.76 


0.04 


3.16 


9.54x10" 


5 


0.148 










Cone(H) 


58 


1.75 


0.40 


3.13 


0.858 




1.04xl0~ 6 


5.97xl0~ 4 








Halo 


102 


2.02 


0.32 


3.73 


0.066 




3.40xl0~ 8 


1.99xl0~ 4 


0.281 






Field 


14 


0.96 


0.29 


3.41 


0.091 




6.60xl0~ 3 


0.065 


0.090 


0.263 




S Mon 


159 


3.05 


1.41 


6.20 
















Spokes cluster 


39 


2.71 


0.72 


6.15 


0.776 














Conc(C) 


74 


2.64 


0.69 


7.54 


0.140 




0.786 








CMD 


Conc(H) 


72 


3.24 


1.67 


8.75 


0.372 




0.824 


0.110 








Halo 


110 


3.80 


1.67 


7.45 


4.05x10" 


4 


0.266 


1.51X10" 2 


0.202 






Field 


33 


3.91 


2.00 


6.50 


2.34x10" 


2 


0.060 


2.82X10 -2 


0.285 


0.922 




Field nm b 


307 


3.66 b 


0.91 b 


10.14 b 


1.04x10" 


4 


0.117 


5.88xl0~ 2 


3.05xl0~ 3 


2.54xl0~ 2 



a Due to the uncertainty in the age of massive PMS stars, the mass limit used here is m = 0.2 - 1.0 Mq. 

b "Field nm" represents the stars in the Field region. They arc in the PMS locus of NGC 2264 CMD, but do not show any appreciable Ha emission or IR 
excess. And therefore their "age" has meaningless because they are not the member of NGC 2264. But to check the reliability of Kolmogorov-Smirnov test, we 
estimate the age of these non-members from the position in CMDs. 
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Table 3. Young Brown Dwarf Candidates 



bung et al. (2008) 


I 


R — I 


V — I 


Ha or X-ray 


2MAsb 


YSO Class 


C3200 


18.290 


1.901 


3.500 


Ha 


06393962+0940442 




C18770 


18.669 


1.977 


3.636 


Ha 


06401412+0934091 




C20186 


18.651 


2.067 


3.946 


Ha 


06401789+0941546 




C21372 


17.508 


2.176 


3.925 




06402107+0935247 


II 


C21459 


19.205 


2.462 




Ha 






C21605 


18.516 


2.394 


4.333 


Ha 






C22265 


18.581 


2.016 


3.540 


Ha 


06402354+0948330 




C23335 


19.023 


2.085 


4.037 


Ha 






C24265 


19.914 


2.365 




Ha 






C24405 


18.554 


2.079 


3.659 


Ha 


06402937+0941035 




C24851 


18.011 


2.184 


3.926 


Ha 


06403063+0954365 


II 


C24999 


17.224 


2.096 


3.840 


Ha 


06403106+0935401 




C25017 


17.525 


2.202 


4.275 


X+Ha 


06403110+0949319 


II 


C28287 


19.907 


2.527 




Ha 


06404043+091 1372 




C28403 


17.739 


2.019 


3.579 


X 


06404078+0934269 




C28559 


18.685 


2.438 


4.326 


x 


06404119+0931287 




C28717 


17.566 


2.050 


3.659 


X+Ha(Sp) 


06404164+0931431 




C28798 


19.147 


2.349 




Ha(Sp) 








17.606 


2.439 


4.452 


X+HaC 


06404494+0938497 




C30192 


17.671 


2.171 


3.931 


Ha(Sp) 


06404610+0947142 




C31370 


20.197 


2.311 




Ha 






C31 392 


17.876 


2.186 


4.169 


X+HafSn") 






C31403 


17.693 


2.129 


3.768 


X 


06405006+0947064 




C31736 


17.549 


2.025 


3.581 


x+Ha 


06405104+0949061 




C31751 


18.152 


1.943 


3.608 


Ha 






C31863 


18.140 


2.013 


4.381 


X 


06405139+0945576 




C32294 


18.829 


2.116 


3.741 


HaC 


06405300+0947377 




C32976 


19.323 


2.407 




X 


06405579+0936061 


II 


C33131 


18.330 


1.973 


3.579 


Ha 


06405641+0934442 




C3321 5 


19.098 


2.323 




Ha 


06405674+09381 01 




C33403 


18.285 


2.481 


4.557 


x+Ha 


06405747+0935147 




C33533 


17.462 


2.102 


3.809 


HaC 


06405798+1002163 


II 


C33590 


17.766 


2.280 


4.122 


X+HaC 






C33624 


17.374 


2.280 


4.256 


X+Ha 


06405834+0937567 


II 


C33627 


20.070 


2.409 




Ha 






C33648 


16.886 


2.363 


4.471 


X 


06405843+0937445 


II 


C33894 


17.675 


2.030 


3.596 


HaC 


06405941+0945516 


II 


C33920 


19.328 


2.207 




Ha 






C33939 


17.641 


2.291 


4.344 


X+HaC 


06405962+0936575 


II 


C34020 


19.809 


2.214 




Ha 






C34101 


18.344 


2.112 


3.813 


HaC 


06410032+0930169 




C34164 


17.430 


2.224 


4.404 


X 






C34674 


18.456 


2.569 




X 


06410258+0936156 




C34695 


17.415 


2.104 


3.774 


X+Ha(Sp) 


06410266+0950329 




C34875 


18.497 


1.995 


3.571 


HaC 


06410350+0935248 




C34965 


18.900 


2.124 


3.820 


X+Ha 
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Table 3 — Continued 



Sung et al. (2008) I R-I V -I Ha or X-ray 2MASS YSO Class 



C34989 


17.661 


2.065 


4.094 


Ha 


06410389+0919352 


II 


C35065 


17.406 


2.296 


4.223 


X+Ha 


06410417+0934572 




C35418 


17.311 


2.227 


4.115 


X+Hq 


06410579+0931072 


II 


C35463 


18.263 


2.343 


4.241 


X+Ha 


06410604+0949232 




C35544 


20.220 


2.301 


4.236 


Ha 






C35632 


18.929 


2.147 


3.939 


X+HaC 


06410696+0935557 




C35763 


18.311 


2.216 


4.994 


X 


06410755+0930004 




C35876 


17.344 


2.401 


4.765 


Ha 


06410812+0918151 


II 


C36048 


17.386 


2.337 


4.963 


X 




II 




1 Q Q9R 






Ha 






C36580 


18.646 


2.508 


4.767 


X+Ha 


06411180+0931123 




C36630 


18.242 


2.071 


4.024 


HaC 


06411208+0945281 


II 


C36702 


17.151 


2.199 


4.069 




06411242+0913247 


II 


C36868 


17.253 


2.302 






06411326+0935473 


II/III 


C36941 


18.357 


2.130 


4.047 


Ha 


06411365+0955370 




C37325 


19.955 


2.548 




Ha 


06411578+0915357 




C37393 


18.666 


2.111 


3.929 


HaC 


06411613+0938166 




C37601 


19.449 


2.178 


3.953 


HaC 






C37602 


18.017 


1.958 


3.523 


Ha 


06411714+0943378 




C37611 


18.672 


1.931 


3.854 


Ha 


06411714+0929048 




C37654 


18.232 


2.257 


4.066 


Ha 


06411736+0931519 




C37667 


19.301 


2.621 




X+Ha 


06411737+0937151 




C37980 


18.770 


2.325 


4.325 


HaC 


06411862+0954074 




C38020 


19.648 


2.315 


4.042 


Ha 






C38517 


19.247 


2.100 


4.313 


Ha 






C38598 


17.402 


2.076 


3.762 


X 


06412104+0932417 




C38760 


18.262 


1.992 


3.806 


Ha 


06412177+1017003 




C39183 


18.760 


2.067 


3.700 


Ha 






C39897 


17.970 


2.047 


3.780 


Ha 


06412593+0930260 




C40179 


18.567 


1.988 


3.528 


Ha 


06412695+0946003 




C42417 


19.487 


2.466 




HaC 






C43063 


18.778 


2.130 


3.976 


Ha 


06413618+1010481 




C44891 


19.875 


2.333 




Ha 







